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^ ; ABSTRACT 

^ ■ The completed 2dF QSO Redshift (2QZ) Survey has been used to search for extreme 

' large-scale cosmological structure ('--^200h~^ Mpc) over the redshift range < z < 2.5. 

. We demonstrate that statistically significant overdensities and underdensities do exist 

' and hence represent the detection of cosmological fluctuations on comoving scales 

CO I that correspond to those presently detected in the cosmic microwave background. 

, However, the fractional overdensities on scales > 100 h~^ Mpc are in the linear or only 

^ ' weakly non-linear regime and do not represent collapsed non-linear structures. We 

\ compare the measurements with the expectation of the ACDM model by measuring 

r-j . the variance of counts in cells and find that, provided the distribution of QSOs on 

C^h' large scales exhibits a mild bias with respect to the distribution of dark matter, the 

Q , observed fluctuations are found to be in good agreement with the model. There is no 

5^ ■ evidence on such scales for any extreme structures that might require, for example, 

I departures from the assumption of Gaussian initial perturbations. Thus the power- 

^ . spectrum derived from the 2QZ Survey appears to provide a complete description 

of the distribution of QSOs. The amount of bias and its redshift dependence that is 
required is consistent with that found from studying the clustering of 2QZ QSOs on 
^lOh"-'^ Mpc scales, and may be adequately described by an approximately redshift- 
' invariant power spectrum with normalisation cts — 1.0 corresponding to a bias at z = 

5^ I of 6 ~ 1.1 rising to 6 ~ 2 at the survey's mean redshift z ~ 1.5. 
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1 INTRODUCTION 

It has long been recognised that active galaxies and QSOs 
may used to probe structure in the universe on the largest 
observable scales, and statistical analysis of the clustering 
of samples of such cosm ological ob.jec t s was first devel- 
oped thirty years ago (e.g.'Webs terllQTdlSeldner fc Feeble^ 
[l978). Even in those early studies it became apparent that, 
although the universe appeared generally homogeneous on 
the largest scal es, individual c lusters of active galaxies could 
be identified llWebsteJ ITosj) . In modern ideas of active 
galaxies and of their link to the process of galaxy forma- 
tion, it is now thought that luminous Q SOs lie at the heart 
of the most massive elliptical galaxi es jKukula et alJl200lt 
iDunlop et al.l2003l:lFiovd et al.l20o3) and that, provided the 
clustering bias expected for such a population is accounted 

* http:/ /www-astro. physics. ox. ac.uk/~Iam 



for JSheth fc TormenligggT) . studying the clustering of active 
galaxies is a means of testing the predictions of hierarchical 
galaxy formation theory both on the largest cosmological 
scales and as a function of cosmic epoch. 

Whether non-uniform distributions of QSOs on large 
scales could in fact be detected with a high degree of statis- 
tical confidence was an open question for many years, how- 
ever, largely because of the limited sample sizes and sys- 
tematic uncertainties and varying selection biases present in 
samples of QSOs. A statistically significant detection of QSO 
clustering on ~10Mpc scales was achieved prim arily from 
the AAT QSO redshift survey JShanks fc BovlJl 9941. With 
the a dvent of the much larger 2dF QSO Redshift Survey 
(2QZ: ICroonTet al. 2001a : Smith et al. 2004; Groom ct al] 
l2004br) not only has an extremely secure detection of QSO 
clustering been achieved, but also the shape of the correla- 
tion function, its bias with respect to the underlying dark 
matter distribution, the cosmic evolution in that bias, and 
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the luminosity dependence of QSO clusterinR have now all 
been measured on scales o ut to about 80 Mpc iCroom et alJ 
l20nihl. I200A l2004al. lioO^ . 

On larger scales, the shape of the power spectrum 
of QSO clustering has been measured and shown to be 
consistent with the expectati ons for a low density CDM 
universe JOutram et al.ll2003l) . The aim of this paper is 
to look at the distribution of QSOs in redshift space, 
rather than Fourier; to investigate on what scales de- 
partures from randomness may be detected; and to test 
whether those departures are consistent with the expec- 
tation of the widely-assumed cosmological model in which 
structures have grown by gravitational collapse of initially 
Gaussian fluctuations in a universe dominated by a com- 
bination of cold dark matter and dark energy. Over the 
years there have been a number of claims of the detec- 
tion of extremely large (50-20 Mpc) groupin gs of AGN and 
QSOs dCram pton ct a UligSTl: Iciowes fc Campusanolll99lt 
IWiUiger et~.i.2002.: .Brand et al.ll2003^ : at the time of these 
claims there was no consensus on the cosmological model 
or on the relationship of the distribution of QSOs to the 
distribution of matter, and it was hoped that the measure- 
ment of such large groupings might constrain models of one 
or both of these. Now, however, there have emerged in the 
literature preference s for one particular cosmological model 
jSoergel et alJ !2003^ and for the bias of QSOs measured by 
the 2dF QS O Redshift Survey (Gr oom et alJ boOlb. 2003 . 
l2004al. |2005|), so it is timely to ask whether the very large 
scale distribution of QSOs is consistent with these recent de- 
velopments. If such large-scale clusters were shown to be col- 
lapsed structures with high values of overdensity this would 
be likely to be in significant disagreement with the cosmolog- 
ical model and measurements of QSO bias. In fact, we shall 
see in this paper that indeed very large groupings of QSOs 
may be identified, possibly on scales as large as 300 Mpc, but 
that the groupings are still in the linear regime of gravita- 
tional collapse and are in fact entirely in accord with current 
expectations. 



2 ANALYSIS OF THE 2QZ SURVEY 
2.1 The 2QZ survey 

The 2QZ survey is described bv [ Groom et al.l ll2001ar) : 
ISmith et all (l2004^ : ICroom et alJ J2004hl) It comprises spec- 
troscopic observations of blue colour-selected QSOs in two 
sky areas totalling 740 deg^: one a strip 5 degrees in dec- 
lination by 75 degrees in right ascension passing through 
the south Galactic pole (the "SGP" region) and the other 
a similar-sized strip on the celestial equator in the region of 
the north Galactic cap (the "NGP" region). In this paper we 
analyse each region separately in order to test for systematic 
errors or differences between the two halves. Spectroscopic 
observations at hj > 18.25 were carried out with the 2dF fa- 
cility at the Anglo- Australian telescope, brighter candidates 
were observed with the 6dF facility at the U.K. Schmidt 
Telescope operated by the Anglo-Australian Observatory 
but are not included in this analysis. The survey contains 
21,181 "Quality 1" QSOs with < 2 < 2.5 and with pho- 
tographic blue magnitude in the range 18.25 < bj < 20.85, 
and this is the sample used in the following analysis. 




Figure 1. The redshift distributions of QSOs in the 2QZ survey, 
binned in redshift intervals of 0.1: (upper) SGP region; (lower) 
NGP region. The smooth curves are the cubic spline fits to the 
observed distributions. 



2.2 The 2QZ selection function and extinction 
correction 

In order to carry out the analysis in this paper it is essen- 
tial to have an accurate measure of the survey's selection 
function. Th i s quest ion has been extensively discussed by 
iGroom et all J2004bil and here we adopt the same empirical 
approach to determining the selection function. 

To do this, we shall assume that the redshift and an- 
gular selection functions are statistically independent and 
hence may be separated - we shall discuss later the validity 
of this assumption. The redshift selection function is de- 
termined simply by fitting a cubic spline to the observed 
redshift distributions in each of the NGP and SGP regions, 
in the range < 2; < 3. The cubic spline was fitted to the 
data binned in redshift intervals of 0.1 with knots at redshift 
intervals of 0.3. Histograms of redshift and the cubic spline 
fits are shown in FigQ The choice of smoothing that results 
from the fitted function may in principle suppress the sig- 
nal arising from real clustering on the largest scales, but by 
averaging over all fields the effects of QSO clustering on the 
redshift distribution should be averaged out, and the result- 
ing function should be a good reflection of the true selection 
function. 

It may be seen that there are some small differences 
between the observed redshift distributions of the NGP and 
SGP regions which may reflect some differences in QSO can- 
didate selection or which may reflect the presence of large- 
scale structure. In this analysis we analyse each region in- 
dependently, using the smoothed redshift distribution from 
each, in order to be able to compare the results from the two 
independent regions and hence to search for any otherwise- 
hidden systematic effects. 

The photometric colour selection used to construct the 
survey becomes inefficient at z > 2.5 JGroom et al.ll2004bl) 
and, as we are concerned about fiuctuations in selection ef- 
ficiency possibly mimicking genuine cosmological structure, 
we have chosen to truncate the analysis at this maximum 
redshift. 

The angular selection function is determined by the dis- 
tribution of photographic plates used for the initial candi- 
date selection, by the magnitude limit and any associated 
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calibration errors, the placement of each spectroscopic 2dF 
field and the configuration of observed targets within each 
field, and the spectroscopic completeness obtained in each 
field. Inspection of the variation in QSO numbers between 
the photographic plates gives no indication of any significant 
calibration errors which would affect the results presented 
here. We shall return to this question below when discussing 
the results. 

To determine the configuration completeness a t any lo - 
cation on the sky, we again follow ICroom et alJ (|2004lJl . 
Each spectroscopic observation was carried out in a 2dF 
field of diameter two degrees. These circular regions were 
overlapped in order to provide contiguous coverage of the 
sky, and the intervals between field centres were chosen to 
optimise the fibre acquisition of targets. Hence we may de- 
fine a set of sectors defined by the boundaries of the 2dF 
fields: a sector may contain data from a single 2dF obser- 
vation or from overlapping 2dF observations, depending on 
the location. The survey comprises 375f such sectors. Within 
each sector, we count the fraction of the colour-selected tar- 
gets that were observed spectroscopically, and define that 
as the configuration completeness at that location. 227 re- 
gions around bright stars were removed from the survey. 
87 percent of the remaining survey area has configuration 
completeness higher than 90 percent. 

However, such a measure does not take into account 
that the efficiency of identifying QSOs may vary between 
2dF fields, depending on observing conditions. Hence an al- 
ternative completeness measure is t he spectroscopic com - 
pleteness, also discussed in detail bv lCroom et al.l ll2004bl) . 
In each sector we determine the fraction of all colour-selected 
targets that have been spectroscopically identified. Hence 
this measure includes the configuration completeness, but 
allows additional variation in spectral quality. 

Although it might seem better to use the spectroscopic 
completeness rather than just the configuration complete- 
ness, there is a concern that this could bias the results, 
because of a combination of two effects. First, it is gener- 
ally easier to identify QSOs than stars from their spectra, 
because of the presence of strong broad emission lines in 
the former. Second, the fraction of targets that actually are 
QSOs varies with sky position, because the surface density 
of stars varies with Galactic latitude and longitude. Hence 
regions of high star density will tend to have lower spec- 
troscopic completeness, even though in fact the efficiency 
of QSO detection may be invariant. This completeness mea- 
sure is also sensitive to variations in the colour selection with 
position: regions with colour selection that is redder, owing 
to photometric calibration uncertainty, would have a higher 
fraction of stellar contaminants, and hence the spectroscopic 
completeness would appear lower, whereas again the QSO 
completeness may be invariant. 

In practice, variations in spectroscopic completeness be- 
tween 2dF fields are unlikely to degrade significantly the 
results presented below, because the angular extent of the 
2dF fields is smaller than the physical scales of interest over 
most of the redshift range, and because neighbouring 2dF 
fields were not typically observed contemporaneously. Hence 
the spectroscopic completeness variations should to some ex- 
tent average out on the scales of interest. For this reason we 
prefer to use the configuration completeness as the baseline 
measure, as this is less susceptible to the other problems 



discussed above, but in the quantitative analysis that fol- 
lows we compare results obtained assuming both complete- 
ness measures, as an indication of the possible magnitude of 
residual systematic uncertainties arising from uncertainties 
in the observational selection. 

Before leaving the discussion of completeness correc- 
tions, we should note that any effects of incorrect complete- 
ness corrections or significant photometric calibration vari- 
ations should be observable as a characteristic signature on 
the maps of fluctuations presented below. A field contain- 
ing more QSOs than its neighbours owing to a complete- 
ness variation would reveal itself as a radial feature in those 
maps. No such features are visible at the level of fluctuations 
detected. 

The final ingredient which should be included as part of 
the selection function is the effect of Galactic extinction. A 
region of high extinction effectively makes the QSO magni- 
tude selection limit brighter, reducing numbers of observed 
QSOs in that region. We correct for the known Galact ic ex- 
tinction using the dust maps of ISchleeel et all (Il998l) with 
resolution of 6 arcmin, an adequate resolution for the detec- 
tion of large-scale fluctuations in QSO numbers. The cata- 
logued values of E{B~V) are converted to extinction assum- 
ing Ab = 4.3I5i5(-B — V) and the effect on QSO numbers is 
calculated assuming a slope for the QSO number-magnitude 
relation at the faint magnitude limit of 0.32. Extinction in 
the NGP region is significantly higher than in the SGP re- 
gion, yet, after application of the extinction correction, no 
statistically significant systematic difference remains in the 
results from the two regions. 

In the analysis presented below the observed fluctua- 
tions in QSO numbers are estimated from a comparison of 
the actual QSO distribution with artificial distributions cre- 
ated from numerical realisations. The completeness and ex- 
tinction corrections are applied to those realisations in order 
to obtain simulated datasets that should mimic the true data 
to high accuracy. 

2.3 The detection of QSO density fluctuations 

A number of proposed methods of detecting clusters of ob- 
jects are described in the literature, but perhaps the most 
straight-forward real-space (or redshift-space) measure in 
terms of understanding its statistics is to smooth the 3D 
distribution of QSOs and to search for statistically signif- 
icant departures from a uniform distribution, taking into 
account sample selection and the shot noise arising from the 
finite sample size. If significant fluctuations on some scale are 
detected then we may further compare the results with the 
expectations of the ACDM cosmological model: this analysis 
is carried out in the next section. In the following analysis we 
assume that redshift-space distortions may be neglected on 
the scales of interest at the amplitudes being investigated. 

As we are interested in trying to detect structures on 
the largest scales that have previously been claimed, we 
smooth the QSO distribution with a spherical top-hat func- 
tion of diameter 100 or 200 h~^ Mpc. The smoothed den- 
sity distribution is evaluated at a sampling of 10 h"'^ Mpc. 
The locations of QSOs are determined assuming the red- 
shift is purely Hubble-flow and assuming a ACDM cosmo- 
logical model with f2j\/ = 0.3, f^A = 0.7. In order to compare 
the observed distribution with the distribution expected in 
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the absence of real density fluctuations, we also compute 
the numbers expected given the survey selection function, 
smoothed and sampled the same as the data. Figs|21&iniplot 
cuts at constant declination through each of the two halves 
of the 2QZ survey for two choices of smoothing scale. The 
plotted scale is the signal-to-noise ratio of the observed vari- 
ations in QSO numbers, estimated as 



^obs 



S ■- 



where riobs is the observed number of QSOs at a given loca- 
tion in the smoothed distribution and Tioxp is the expected 
number given the selection function and in the absence of 
any cosmological structure. Only regions with riexp > 10 are 
plotted. The statistic S may only be interpreted as hav- 
ing a normal distribution in the limit of large values of 
riexp, but in fact given the typical value of ncxp ~ 100 this 
is a reasonable approximation, and leads to the observa- 
tion that when smoothed on this scale, there exist statis- 
tically significant deviations from the null hypothesis that 
QSOs trace a randomly sampled uniform distribution. The 
regions that appear significantly under- or over-dense are 
resolved on these maps, and it appears that significant ex- 
cesses may be traced out to substantially larger scales than 
the 200 h~^ Mpc smoothing scale. Although the top-hat fil- 
ter has undesirable properties when transformed to Fourier 
space, in the analysis presented here one can be confident 
that distant points in the smoothed distribution have un- 
correlated shot noise, and the amplitude of fiuctuations ob- 
served in the smoothed maps may be more readily compared 
with the analysis of the variance of counts in cells that is 
presented in the next section. 

We should perhaps be concerned that errors or uncer- 
tainties in the 2QZ selection function may result in spurious 
creation of fluctuations in Figs|5|&|3| The strongest argu- 
ment against this possibility are that the selection function 
and extinction correction have been entirely separated into 
radial and angular components, and we would expect errors 
in either component to be manifest as purely radial or an- 
gular signals. In Figs|5|&|21 we also show the boundaries de- 
fined by the right ascension limits of the UKST photographic 
plates used to construct the survey. There is no obvious ten- 
dency for there to be radial overdensities aligned with these 
boundaries: although the alignment of a minority of the fea- 
tures appears radial, this is not generally the case, and those 
features that could be radial in nature may be seen to not fill 
an entire UKST photographic plate and to extend across the 
plate boundaries. Overall, we conclude that there are no ob- 
vious radial signals associated with the photographic plates. 
At a greater level of subtlety, there may be variations in the 
selection function which are not well represented by the as- 
sumption of separable components (i.e. angular-dependent 
variations in the redshift selection function). However, if 
such effects are present, the existence of under- and over- 
dense regions apparently at randomly placed redshifts, and 
not particularly restricted to individual photographic plates, 
implies that such an effect is not dominating the observed 
distributions. 

Figs^l&El show the same data but now presented as 
plots of fractional overdensity in QSO number. Because 
shot noise cannot be neglected, for this exercise we adopt 
a Bayesian estimator of the fractional overdensity at each 



location on the map, which maximises the posterior proba- 
bility 

P (5i"-obs, Wcxp, O-^) OC P {SW^) P ("-obsi^, ncxp) 

where 5 is the fractional overdensity and a'^ is the cosmologi- 
cal variance in 5 on the scale being considered. For Gaussian 
fluctuations 
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The Bayesian estimate has the property that in the presence 
of random noise the expectation value of a set of points with 
the same estimated value is the true value, and in that sense 
it provides the "most likely" value of fluctuation at any point 
on the map. A maximum likelihood estimate would have the 
property that the expectation value of a set of points with 
the same true value is that value. Taking the maximum like- 
lihood estimate in the presence of random noise leads to the 
statistical distribution of fluctuations being broader than the 
true distribution, whereas the Bayesian estimator leads to a 
distribution that is narrower than the true distribution. On 
the maps presented here, regions of low signal-to-noise have 
Bayesian estimates of overdensity that are biased towards 
zero (if a maximum likelihood estimate had been applied in- 
stead, regions of low signal-to-noise would appear biased to 
values more extreme than reality). This effect is noticeable 
in Figs|l]&|^at low redshift, where the increasing shot noise 
leads to some dilution of the measured signal. The Bayesian 
estimates were created assuming the ACDM model of Sec- 
tion 2.4 with as — 1.0. 

It is immediately seen that the detected fluctuations are 
all in the linear or weakly non-linear regime, with —0.31 < 
S < 0.44 on scales of 100 h"^ Mpc and -0.13 < S < 0.15 on 
scales of 200 h~^ Mpc. Hence the chief result of this section 
is that large-scale fluctuations in QSO numbers exist, as sus- 
pected by previous authors, but unless the bias is much less 
than unity (the issue of QSO bias is discussed later) the mass 
density fluctuations are not of sufficiently large amplitude to 
represent collapsed structures. Hence the best matching hy- 
pothesis is that the fiuctuations are simply those that are 
expected in the linear regime of the growth of cosmological 
structure, and that these are fully represent ed statistically 
by th e power spectrum of the 2QZ survey llOutram et alJ 
l2003r . This hypothesis is tested further in the next section. 

2.4 A counts-in-cells analysis 

Having demonstrated that statistically-significant fiuctua- 
tions do exist on large scales, we now should test how these 
fiuctuations compare with the predictions of the cosmologi- 
cal model assuming that on large scales structures are fully 
represented by the power spectrum. To do this we carry 
out a counts-in-cells analysis, measuring the sample vari- 
ance when the data are sampled in independent cubes of a 
specified size, and comparing that sample variance with the 
variance expected given a particular choice of power spec- 
trum and normalisation. 

The measured cell-to-cell variance has contributions 
from two sources: cosmological fiuctuations and shot noise. 
Furthermore, because the QSO density varies over the sur- 
vey, the weight given to each cell should be varied to op- 
timise the signal-to-noise of the overall sample variance. In 
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Figure 2. Cuts at constant declination tlirougli tlic sniootlicd tlircc-dinicnsional maps of fluctuations in the NGP (upper) and SGP 
(lower) regions. The x-axis corresponds crudely to the right ascension direction, the y-axis to the redshift direction. The smoothing 
function is a spherical top hat of diameter 100 h^^ Mpc as indicated on the lower left of each map. The grey scale is in units of signal-to- 
noise ratio S, as defined in the text, and is the same on both plots. The greyseale bar indicates the maximum and mininnun values found 
in each region. Also shown are the boundaries defined by the right ascension limits of the UKST photographic plates used to construct 
the survey. 
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Figure 3. As Fig|^ but with a splierical top liat smoothing function of diameter 200 h Mpc as indicated on the lower left of each 
map. 
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Figure 4. Cuts at constant declination tliiougli tlic Baycsian maps of smoothed fluctuations in the NGP (upper) and SGP (lower) 
regions. The smoothing function is a spherical top hat of diameter 100 fi^^ Mpc as indicated on the lower left of each map. The grey 
scale is in units of fractional overdensity and is the same on both plots. The bar also indicates the maximum and minimum values found 
in each region. 
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Figure 5. As FigBlbut with a splierical top hat smoothing function of diameter 200 h~^ Mpc. 
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the limit where QSO numbers are sufficiently large that the 
shot-noise Poisson distribution tends to a normal distribu- 
tion, the maximum likelihood estimate of variance of over- 
density S is 

2 y ] (^obs ^cxp) y ] '^cxp 

^ ~ 2 ■ 

As here rioxp >> 1 we adopt this statistic. However, as com- 
parison is made with values of this statistic measured in 
simulated data (see below) we do not need to assume that 
the central limit theorem holds. The statistic is calculated 
from the data for cell sizes of 80, 160, 320 & 640 h"^ Mpc. 
The use of cubic cells means that, for a given cell size, each 
cell has shot noise that is statistically independent from its 
neighbour (although the cosmological signal does not have 
this property). The variances we expect then are slightly dif- 
ferent from the variances that would be calculated using a 
more standard assumption of spherical cells. 

In order to calculate the expected variance for a given 
cosmological model, we could in principle carry out the 
standard integration over the power spectrum (e.g. iPeacocld 
ll99S) assuming an appropriate cell size and shape. There are 
a number of reasons not to do this, however. First, the sur- 
veyed volume does not fill a geometrically simple region, and 
in fact on the largest scales probed here the survey becomes 
two-dimensional at low redshifts: the survey dimension in 
the declination direction is smaller than the largest cell size. 
This causes the expected variance to have a modified depen- 
dence on cell size. Second, the fiuctuations in cells of differ- 
ing sizes are correlated, because the same data are used in 
the analysis of different cell sizes. The degree of correlation 
depends on the selection function, the power spectrum, the 
number of QSOs and the choice of statistic used to measure 
the variance. 

Hence the approach adopted here is to create mock QSO 
catalogues which have the same selection function and shot 
noise as the actual data, and in which the expected level 
of cosmological fluctuations are present, assuming a partic- 
ular cosmological model. By averaging the values obtained 
from a large number of such simulations this approach au- 
tomatically includes the effects of cosmic variance in the 
distribution of the statistic s^. The simulations are gener- 
ated by creating random density fields on a grid of sampling 
10 h~^ Mpc and size much larger than each region of the 2QZ 
survey, whose Fourier components are drawn randomly from 
a normal distribution with variance specified by the power 
spectrum. The density fields are then multiplied by the 
2QZ selection function and extinction corrections. These are 
Poisson-sampled to create simulated datasets whose proper- 
ties should mimic in every way the properties of the actual 
data. This process is repeated for the NGP and SGP regions 
of the survey. 

The shape of po wer spectrum adopted is that of 
lEfstathiou et al.l ^9^. We assume a ACDM cosmology 
with f^A = 0.7, Qm = 0.3, Ho = 70kms"^Mpc"\ To set 
the normalisation we need to consider the bias of the power 
spectrum, but dealing with this is complicated by the possi- 
bility that departures of the Poisson shot noise from a nor- 
mal distribution could lead to an offset between the values 
of the statistic measured for the data and for the simu- 
lations. Hence we have tested a range of values of the nor- 
malisation parameter erg (the rms fluctuation in spheres of 



diameter 8h~^ Mpc), and choose the value which produces 
the best fit of the statistic to the data. The dark mat- 
ter power spectrum also evolves in amplitude with redshift: 
what should we assume for the power spectrum of QSO num- 
ber fiuctuations? In fact, we expect QSOs to exhibit larger 
fiuctuations than the dark matter, because QSOs exist in 
the m ost massive individual galaxies known jPunlop et al.l 
l2003h . and hence the amplitude of QS O fluctua tions will be 
biased, as described by, e.g. Shet h fc TormenI (_1999). Be- 
cause massive galaxies are increasingly rare at higher red- 
shifts, the bias increases with redshift. Hence the expected 
amplitude of QSO fluctuations should have a redshift depen- 
dence which is dependent on the mass of QSO host galaxies. 
Rather than assuming a theoretical value for the evolution of 
the bias, we ca n measure it from the 2QZ survey on smaller 
physical scales dCroom et alJl2001bll2004al . l2005l. Loaring et 
al. in prep.) and assume that the cosmological evolution of 
the bias is the same on large scales as on 10 Mpc scales. To 
a good approximation, in a ACDM cosmology, the lOMpc- 
scale clustering has an amplitude which is constant with 
redshift over the range 0.5 < z < 2, so here we shall as- 
sume that the power-spectrum amplitude is invariant with 
redshift. The uncertainties in the variance of the cell counts 
are large, so it should make little difference if in reality there 
is some departure from the "no-evolution" assumption (in 
fact, it would be more correct to say that there is strong 
evolution of the bias which cancels the evolution in the dark 
matter fluctuations). 

Fig-ini shows the measured variance as a function of 
cell size for the two halves of the survey. It may be seen 
that the two halves are in good agreement. Values are 
shown for both methods of completeness correction: config- 
uration (filled symbols) and spectroscopic (open symbols). 
Also shown on the figure is the expected variance from the 
mean of 100 ACDM simulations for a QSO power spectrum 
normalisation of ag — 1.0, this being the best fitting redshift- 
independent power spectrum normalisation. Extrapolated 
to z = this would correspond to a bias value b ~ 1.1, 
but at the typical redshift of the 2QZ sample the bias 
is substantially higher, consistent with the findings based 
on measuring the correlation f unction on smaller scales 
JCroom et alJl200ll2004all2005ft . The error bars shown re- 
fiect the 68 percent range of values in cell variance obtained 
from the simulations and are therefore good estimates of the 
random errors, taking into account both shot noise and cos- 
mic variance. The total of the fit to both regions is 1.8 
with 7 degrees of freedom, assuming the selection function 
determined from the configuration completeness. The covari- 
ance between data points in Fig.|S|is low but is included in 
the estimation. The measured variance increases and the 
best-fitting value rises slightly to 5.3 if the spectroscopic 
completeness correction is adopted (open symbols in Fig.jSJ. 
It may be seen that on the largest scales the two complete- 
ness correction methods reveal systematic uncertainties in 
the variance that are comparable to the size of the random 
errors (shot noise and cosmic variance). 

If we assume as — (i.e. that there are no cosmological 
fiuctuations) rises to 61, indicating that the hypothesis 
that the observed fiuctuations are merely due to shot noise 
may be rejected at a significance level of 3 x lO"'^". 

We note at this stage that we have not attempted to 
compare the data with any competing models, but rather 
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Figure 6. The variance of counts in cells for the NGP and SGP halves of the survey. Filled points show the measured results after 
correcting for configuration completeness. Open points show the results after correcting also for spectroscopic completeness. Also shown 
are mean values and 68 percent confidence intervals derived from the numerical simulations for each region and for each of the two 
completeness corrections. A value of a-g = 1.0, which best fits the filled points, has been adopted. 



simply to test whether the data are consistent with the 
model favoured by independent cosmological measurements: 
the approach adopted here may be viewed as being frequen- 
tist in nature rather than Bayesian. However, we might won- 
der whether the shapes and topology of the detected struc- 
tures are as expected within this model, or whether there 
is any evidence for significant differences. Statistical tests of 
topology are beyond the scope of this paper, but by compar- 
ing the structures detected on the two scales shown in Figs 2- 
5, we can see qualitatively at least that there is no evidence 
for, for example, extremely filamentary structure being re- 
sponsible for the observed fluctuations on the larger scale. A 
qualitative comparison of the appearance of Figs 4 & 5 with 
the ACDM simulated maps (not shown here) also reveals no 
significant difi^erences. 



3 DISCUSSION AND CONCLUSIONS 

Fluctuations in QSO space density on scales ~ 200 h~^ Mpc 
may be seen directly in the QSO distribution (Figs l2l5Tl . 
It is likely that these fluctuations may be traced out to 
larger scales, and inspection of the maps and analysis of the 
variance of the QSO density field indicates the detection 
of fluctuations on scales possibly as large as 300 h~^ Mpc. 



The fluctuations are in good agreement with those expected 
in a ACDM cosmology with WMAP parameters if we as- 
sume that QSO clustering is a biased tracer of dark matter 
fluctuations, with a bias value that is approximately red- 
shift independent and which is close to the value inferred 
from QSO clustering on 10 Mpc scales. Since it appears that 
6 > 1, the fluctuations in mass associated with the observed 
QSO fluctuations are certainly in the linear or weakly non- 
linear regime of gravitational collapse, and hence they do not 
represent collapsed non-linear structures, but are merely a 
reflection of the large-scale fluctuations expected given the 
ACDM power spectrum. Figs l2l5l represent direct detection 
of structure in the distribution of discrete objects on comov- 
ing scales that correspond to the scales measured in the cos- 
mic microwave background by WMAP ( SpcrEcl ot al. 200^. 

Consideration of two methods of correcting the mea- 
surements for observational selection shows that on the 
largest scales there may remain systematic variations in sur- 
vey uniformity at levels that are comparable to the measured 
fluctuations. For this reason we believe that the 2QZ survey 
should not be used to attempt to infer more precise values 
for cosmological parameters such as Q.m until those system- 
atic non-uniformities can be better corrected for. 

One of the most signiflcant features of the results pre- 
sented here is that there is no evidence for any non-Gaussian 
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initial conditions, which might have been implied had any 
of the > 50 Mpc structures previously claimed for QSO dis- 
tributions turned out to represent collapsed structures. As 
the counts in cells are completely in accord with the simula- 
tions generated from a ACDM power spectrum, we conclude 
that, as far as this test is concerned, the power spectrum re- 
mains a complete description of the large-scale distribution 
of QSOs. 

We have not in this paper attempted to measure the 
topology of the oljserved fluctuations. This would be an 
interesting exercise in order to test whether the observed 
structures are consistent in every way with the expectations 
of linearly collapsing structures forming from Gaussian ini- 
tial conditions. 
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